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1. Introduction 

It is presently widely accepted that quantitative understanding of confinement of quarks 
and gluons remains as the major challenge in QCD. In this connection, a linearly 
increasing quark-antiquark pair static potential provides the simplest criterion for 
confinement, although unfortunately there is up to now no known way to analytically 
derive the confining potential from first principles. In this context, it may be recalled 
that phenomenological models have been of considerable importance in order to provide 
strong insight into the physics of confinement, and can be considered as effective theories 
of QCD. One of these, that is the dual superconductivity picture of QCD vacuum P, 
has probably enjoyed the greatest popularity. The key ingredient in this model is the 
condensation of topological defects originated from quantum fluctuations (monopoles). 
As a consequence, the color electric flux linking quarks is squeezed into “strings”, and 
the nonvanishing string tension represents the proportionality constant in the linear, 
quark confining, potential. We further note that recently an interesting approach to 
this problem has been proposed [ 2 ], which includes the contribution of all topologically 
nontrivial sectors of a gauge theory. Mention should be made, at this point, to lattice 
calculations which clearly show the formation of tubes of gluonic fields connecting 
colored charges p. In agreement with lattice results, loop-loop correlations have been 
recently computed analytically in extended stochastic vacuum models |lj. 

With these ideas in mind, in a previous paper j3], we have studied a simple effective 
theory where confining potentials are obtained in the presence of nontrivial constant 
expectation values for the gauge field strength coupled to a scalar (axion) field 
via the interaction term 

= ( 1 ) 

O 

In particular, we have observed that in the case of a constant electric field strength 
expectation value the static potential remains Coulombic, while in the case of a constant 
magnetic field strength expectation value the potential energy is the sum of a Yukawa 
and a linear potential, leading to the confinement of static charges. More interestingly, 
we remark that the magnetic character of the field strength expectation value needed 
to obtain confinement is in agreement with the current chromo-magnetic picture of 
the QCD vacuum [B]. Another feature of this model is that it restores the rotational 
symmetry (in the potential), despite of the external fields break this symmetry. We 
further observe that similar results have been obtained in the context of the dual 
Ginzburg-Landau theory [7], as well as, for a theory of antisymmetric tensor fields 
that results from the condensation of topological defects as a consequence of the Julia- 
Toulouse mechanism jH]. Accordingly, from a phenomenological point of view, we 
have established an equivalence between different models describing the same physical 
phenomena. This allows us to obtain more information about a theory than is possible 
by considering a single description. 

By following this line of reasoning, it is natural to extend the previous analysis to the 
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case where a scalar field 0 is coupled to a non-Abelian gauge field, that is, 

= . ( 2 ) 

O 

where the trace is taken over color indices. While an Abelian model with interaction of 
the type (CD is no more than a useful “laboratory” , it may be worth to recall that in 
QCD coupling of the form m is instrumental to build up the Peccei -Quinn mechanism 
[T^ and solve the the strong CP problem p cni El. In this case, the scalar field 
describes the axion, i.e. the Nambu-Goldstone boson of a new broken 17(1) symmetry 
of the quark and Higgs sector. 

The purpose of this work is to extend the Abelian calculations in P to the non-Abelian 
case and find the corresponding static potential. Our calculations are done within 
the framework of the gauge invariant/path-dependent variables formalism, providing an 
effective tool for a better understanding of effective non-Abelian theories. One important 
advantage of this approach is that provides a physically-based alternative to the usual 
Wilson loop approach, where in the latter the usual qualitative picture of confinement 
in terms of an electric flux tube linking quarks emerges naturally. As we shall see, in the 
case of a constant chromo-electric field strength expectation value the static potential 
remains Coulombic. On the other hand, in the case of a constant chromo-magnetic field 
strength expectation value the potential energy is the sum of a Coulombic and a linear 
potential, that is, the confinement between static charges is obtained. As a result, the 
new coupling displays a marked departure of a qualitative nature from the results of 
Ref. P at large distances. 

It is interesting to observe that the dual superconductivity picture of QCD, compared 
with the model proposed here, involves the condensation of topological defects originated 
from quantum fluctuations. Thus one is led to the conclusion that the phenomenological 
model proposed here incorporates automatically the contribution of the condensate 
of topological defects to the vacuum of the model or, alternatively, the nontrivial 
topological sectors as in p. Another way of obtaining the above conclusion is by 
invoking the bosonization technique in (l-l-l)-dimensions. Indeed, it is a well known fact 
that, for instance, in the Schwinger model na. the bosonized version (effective theory) 
contains quantum corrections at the classical level. In the same way, we can interpret 
the model proposed here as an effective theory which contains quantum effects at the 
classical level. Thus, one obtains a similarity between the tree level mechanism that 
leads to confinement here and the nonperturbative mechanism which gives confinement 
in QCD. Accordingly, the above interrelations are interesting from the point of view 
of providing unifications among diverse models as well as exploiting the equivalence in 
explicit calculations, as we are going to show. 

2. Interaction Energy 

As we discussed in the introduction, our immediate objective is to calculate explicitly 
the interaction energy between static pointlike sources, for a model containing the term 
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along the lines of Refs. Emni. To this end we will compute the expectation value 
(H)^ of the Hamiltonian operator H in the physical state |<h) describing the sources. 
The non-Abelian gauge theory we are considering is dehned by the following generating 
functional in four-dimensional spacetime: 

^ = J TtfiVAexp i^—i J d'^xC^ , (3) 

with 

C = -^TrF^,F^" + TrF^.F,. + , (4) 

where is the mass for the axion field fi. Here, (x) = Af (x) Ta, where 
Ta is a Hermitian representation of the semi-simple and compact gauge group; and 
= d^Af — duAf^ + g ff^A^^Af, with the structure constants of the group. As in j3] 
we restrict ourselves to static scalar helds, a consequence of this is that one may replace 
d'^cj) = —V^0. It also implies that, after performing the integration over in Z, the 
effective Lagrangian density is given by 

£ = -1 TVF^„F'"' + . (6) 

Furthermore, as was explained in jT3], the expression © can be rewritten as 

C = -^Trf^^R" + ^e/^-Z^Tr (F^,) ^ ' ^ 6 ) 

where, represents the constant classical background (which is a solution of 

the classical equations of motion). Here, describes a small fluctuation around 
the background, we also mention that the above Lagrangian arose after using 
£^'^“/^Tr (F^j,) (fFag) = 0 (which holds for a pure chromo-electric or a pure chromo- 
magnetic background). 

By introducing the notation and {Ppo) = expression (jHI) 

then becomes 

C = -^Trf^^R^ + ^Trv“/3f«^ ^ ^ • (7) 

4 16 ^ V"' — m\ 

At this stage we note that m has the same form as the corresponding Abelian effective 
Lagrangian density. This common feature is our main motivation to study the effect of 
the non-Abelian coupling on the interaction energy. 


2.1. Chromo-magnetic case. 


We now proceed to obtain the interaction energy in the v“ 7 ^ 0 and v*-^ = 0 case 
(referred to as the chromo-magnetic one in what follows), by computing the expectation 
value of the Hamiltonian in the physical state |<F). The Lagrangian ((7j) then becomes 

£ = + ^Tlv“f„=5-!_ 

4 4 — m\ 


- TrAoJ° 


( 8 ) 
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where J° is an external cnrrent, (/x, z/ = 0,1, 2, 3) and {i, k = 1,2,3). Once this is done, 
the canonical qnantization in the manner of Dirac yields the following resnlts. The 
canonical momenta are 


= 0 , 

(9) 

-pra _ e^ab p j 

(10) 

III 

(11) 

/ /32 1 \ 

TACib _ I i ^ ^,a ^,6 \ 

A, -(^.5 A + 2''“v2-mT)7 ' 

(12) 


Since D is a nonsingular, there exists the inverse D ^ and from Eq. m we obtain 


Ef = 


det D 


det D - ^ 




2 




(13) 


The corresponding canonical Hamiltonian is thus 


He = I dYTi 


nODAo), + -n*ni + -B‘B, 




V 2 -M 2 


(14) 


where 


—and B* is the chromo-magnetic held. By applying Dirac 
quantization procedure for constrained systems, and removing non-physical variables 
by imposing an appropriate gauge condition ||, we can compute the interaction energy 
between pointlike sources in the model under consideration. A fermion is localized at 0 
and an antifermion at y. From our above discussion, we see that {H)^ reads 

1 


{H)^ = ($1 / d^xDi 


1 • <52 . 

-Trffn,- - —Trv*n,- 


V 2 - M 2 


Trv*n,- 


Id)) .(15) 


Now we recall that the physical state can be written as |15j . 

\Y dz^Ai (z)'^ (0) \0) . (16) 

The line integral is along a spacelike path on a hxed time slice, P is the path¬ 
ordering prescription and |0) is the physical vacuum state. As in jin], we again restrict 
our attention to the weak coupling limit. From this and the foregoing Hamiltonian 
discussion, we then get 


{//>, = {//>„ + V, + 14 , (17) 

where {H)^ = (0| M |0), and the Vi and V 2 terms are given by: 

Hi = ^(®l/d^xTrn^n, Id)), (18) 

II In our gauge-invariant, path-dependent formalism, gauge fixing procedure is equivalent to the choice 
of a particular path US], e.g. a spacelike, straight, path = P X {x — pf, on a fixed time slice. 
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and 


V2 = -^{^\J d^xWUi 


■Trv*ni|d)) 


(19) 


V 2 - M 2 

One immediately sees that the Vi term is identical to the energy for the Yang-Mills 
theory. Notwithstanding, in order to put our discussion into context it is useful to 
summarize the relevant aspects of the analysis described previously ng. Thus, we then 
get an Abelian part (proportional to Cp) and a non-Abelian part (proportional to the 
combination CfCa)- As we have noted before, the Abelian part takes the form 

= l-g^TiiT^^Ta) r dC r dz, 6 (z - z') , (20) 

2 JO JO 

remembering that the integrals over and z^ are zero except on the contour of 
integration. Writing the group factor TrT“Ta = Cp, the expression (P|) is given by 

(i) = ~ . ( 21 ) 

where |y| = L. Next, the non-Abelian part may be written as 

=Tt J d^x{0\l%\0) , (22) 

where, 

1“ = jV,” r‘ r dz^ f dXAt (z) z‘d (x - Az). (23) 

JO JO 

It should be noted that, by using spherical coordinates, expression (121 reduces to 

1 


jai _ 2 fa rpb 

^ — y JbC TTl i2 

Z X 


dz'‘At{z)j2yL{o',^')yi,n{i),^) 


Im 


Putting this back into Ea. (|22|l . we obtain 


{L) = -CaCp 


d)/.’ 


dY / dz'^Dij {z, z') 


(24) 


(25) 


Here Dij{z,z') stands for the propagator, wich is diagonal in color and taken in an 
arbitrary gauge. Following our earlier discussion, we choose the Feynman gauge. As a 
consequence, expression (121 then becomes 

(L) = - 9 *log (AL) , (26) 

where A is a cutoff. Then, the Vi term takes the form 

V, = -g^Cp^(^l + ^CAlog{AL)j . (27) 

It is important to realize that our calculation was based only on the Hamiltonian and 
on the geometrical requirement that the fermion-antifermion state be invariant under 
gauge transformations. From (El we see that the term of order g'^ is just the Coulomb 
energy due to the color charges of the quarks. The correction term of order g‘^ represents 
an increase of the energy due to the vacuum fluctuations of the gauge fields. 
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The task is now to evaluate the V 2 term, which is given by 

-Trv^nj|<h) . 


14 = -T 


(28) 


Once again, from our above Hamiltonian structure we have an Abelian contribution and 
a non-Abelian contribution, in other words. 


vr-' = / dV £ d-J'S (X - z') —^ £ 


dz^6 (x — z) , 
(29) 


and 




(Non—Ab.) _ 


•) = ^Tr ( ) UcAX'vf [ dz^ £ dz'^Dik (z, z') £ du^ 


as before, Dik{z,z') represents the propagator. Here, G is the Green function 

1 g-M|z'-z| 


G (u, v) = 


dvr |z' 


This Green function is, in momentum space, 
1 r d^k 


(31) 


(32) 


dvr |u —v| J (27r)^ 

By means of Eq. (IH^ and remembering that the integrals over z* and z'^ are zero except 


on the contour of integration, the term (i2nD reduces to the linearly increasing potential 
|2j, that is. 




(33) 


We now proceed to calculate the term. As before, we will use the Green 

function in momentum space to handle the integral in Eq. dsni). Following our earlier 
procedure j3], Ea. dHIH) is further rewritten as 


yiNon-Ab.) ^ X 


f 1 + dz^ £ dz"^ |z| Dik (z , z') 


(34) 


Now, we move on to compute the integral (IHH) . As in the previous calculation, we 
choose Dik{z,z') in the Feynman gauge. Thus the jg 

yyn-Ai.) = !t£ln( TV) r'^gyviL log ii + , (35) 

after substracting the self-energy terms. 

From and dsg) we then get 


2 z ,2 r 


W = 


8 


^Tr ( vV,) + /Tr ( T^T'^) Ubcfl,v^'v\ 
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By putting together Eqs. (EZD and (ESI), we obtain for the total interquark potential 

,2 


V= - 


9^Cf 

AttL 


l + ^C'^log(AL) 

TT 


+ 


^ ^ ^ ^Tr(vVi) + Tr ( 


fpbcfqd 




L\og 




(37) 


It must be observed that the rotational symmetry is restored in the resulting form of 
the potential, although the external helds break the isotropy of the problem in a manifest 
way. It should be remarked that this feature is also shared by the corresponding Abelian 
interaction energy [S]. As we have noted before this improves the situation as compared 
to the ’’spaghetti vacuum” model [S] where rotational symmetry seem to be very difficult 
to restore. 

Now we recall the calculation reported in [2j by taking into account topological nontrivial 
sectors in 17(1) gauge theory 

We immediately see that the result E7|) is GXcLctly tliG oiiG obtciinGci by R,Gf. [2j* It is 
interesting to notice that even if is an effective model, extending the Abelian one 
discussed in j^, it is able to reproduce the correct form of the Cornell potential 
As such it deserves some further investigation. 


2.2. Chromo-electric case. 

Now we focus on the case = 0 and 7 ^ 0 (referred to as the chromo-electric one in 
what follows). The corresponding Lagrangian density reads 

C = -tTrf^^P" + ^Trv*^f,,—- TrAoJ° , (39) 

4 16 ■' — m\ 

(/i, u = 0,1, 2, 3) and {i,j, k,l = 1, 2, 3). 

Here again, the quantization is carried out using Dirac’s procedure. We can thus write 
the canonical momenta = —/“°^, which results in the usual primary constraint 
= 0 and Dehning the electric and magnetic helds by and 

Qak _ jaij^ respectively, the canonical Hamiltonian is thus 


He = I d^x 


1 • 1 ■ 
-TrE*Ei + -TrB'Bi - — 
2 2 16 


£pm£fc/nTrv*^B^ 


V2 


-Trv'^^B” 


m. 




(40) 


It is straightforward to see that the constrained structure for the gauge held is identical 
to the usual Yang-Mills theory. However, in order to put the discussion into the context 
of this paper, it is convenient to mention the relevant aspects of the analysis described 
previously ng. Therefore, we pass now to the calculation of the interaction energy. 












Confinement effects from interacting chromo-magnetic and axion fields 


9 


As done above, our objective is now to calculate the expectation value of the Hamiltonian 
in the physical state |<h). In other words, 


= ($1 - d^xTiE^E, 1 $) 


(41) 


' ' 2 ^ 

Taking into account the above Hamiltonian structure, the interaction takes the form 
{H)^ = {H)o + V^, (42) 

where {H)q = (0 | if | 0). Accordingly, the potential reads 




(43) 


3. Final Remarks 


is chromo-electric-like no unexpected features 


In summary, we have considered the conhnement versus screening issue for a non- 
Abelian theory with a coupling in the case when there are nontrivial 

constant expectation values for the gauge held strength The constant gauge held 
conhguration is a solution of the classical equations of motion. 

It was shown that in the case when 
are found. Indeed, the resulting static potential remains Coulombic. More interestingly, 
it was shown that when is chromo-magnetic-like the potential between static 

charges displays a Coulomb piece plus a linear conhning piece. An analogous situation 
in the Abelian case may be recalled [H] . Also, a common feature of these models (Abelian 
and non-Abelian) is that the rotational symmetry is restored in the resulting interaction 
energy. 

We recall that the ehective action m is the expansion of Q up to second order in 
the huctuation held thus one could wonder about the ehect of including higher 
order terms. Even if an explicit calculation is well beyond the purpose of this note, 
some general comments can be done. We know that higher order quantum ehects 
renormalize coupling constants by making them scale dependent. It is easy to see that 
the logarithmic correction to the Coulombic term in (EZD is nothing but the hrst order 
expansion of the ehective, running, coupling constant 


alff 


AL) = 


g'^TT 


l-^log(AL)’ 


(44) 


gljrj- exhibits the expected asymptotically free behavior at short distance 
characterizing the non-Abelian character of the strong interaction. This result marks 
a clear diherence from the Abelian case, where the electric charge increases at short 
distance. Similar renormalization ehects are expected for the string tension a, as well. 
However, the log{l -|- ■^) factor in the second term in (pTTj) is L-independent. We 
are conhdent that higher order corrections will preserve this feature leading to a still 
conhning static potential of the form 


nr) = ^ 


4L 


(45) 
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An explicit check of (gSI) will be addressed in a future work. 

We conclude by noting that our result agrees with the monopole plasma mechanism 
HI. a However, although both approaches lead to conhnement, the above analysis 
reveals that the mechanism of obtaining a linear potential is quite different. As 
already mentioned, in this work we have exploited the similarity between the tree level 
mechanism that leads to conhnement here and the nonperturbative mechanism (caused 
by monopoles) which gives conhnement in QCD. 
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